Kondo effect in a quantum dot coupled to ferromagnetic leads and side-coupled to a 

nonmagnetic reservoir 
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Equilibrium transport properties of a single-level quantum dot tunnel-coupled to ferromagnetic 
leads and exchange-coupled to a side nonmagnetic reservoir are analyzed theoretically in the Kondo 
regime. The equilibrium spectral functions and conductance through the dot are calculated using 
the numerical renormalization group (NRG) method. It is shown that in the antiparallel magnetic 
configuration, the system undergoes a quantum phase transition with increasing exchange coupling 
J, where the conductance drops from its maximum value to zero. In the parallel configuration, on 
the other hand, the conductance is generally suppressed due to an effective spin splitting of the dot 
level caused by the presence of ferromagnetic leads, irrespective of the strength of exchange constant. 
However, for J ranging from J = up to the corresponding critical value, the Kondo effect and 
quantum critical behavior can be restored by applying properly tuned compensating magnetic field. 

PACS numbers: 72.25.-b, 73.63.Kv, 73.23.-b, 73.43.Nq, 85.75.-d 



I. INTRODUCTION 

Transport through a model single-level quantum dot 
captures many interesting and important features of 
transport phenomena in real quantum dots. One of such 
phenomena, which has been of great interest in the last 
decade, is the Kondo effect^— When the dot is occupied 
by a single electron, virtual transitions between the dot 
and electron reservoirs (external leads) cause spin fluc- 
tuations in the dot. As a result, the dot's spin becomes 
screened by electrons of the reservoirs, which results in 
the formation of a non-local spin singlet ground state of 
the system. Furthermore, a resonance in the density of 
states appears at the Fermi level, which gives rise to en- 
hanced transmission through the dot. In experiments, 
this leads to the well-known zero-bias anomaly, i.e. a 
peak at zero bias in the differential conductance*^ 

When the reservoirs are ferromagnetic, the effective 
exchange field generated by the electrodes may sup- 
press the Kondo anomaly^ - — More specifically, when 
the dot described by an asymmetric Anderson model is 
symmetrically coupled to ferromagnetic leads, then the 
Kondo effect becomes suppressed in the parallel configu- 
ration, while in the antiparallel configuration the Kondo 
anomaly survives. However, the Kondo effect in the 
parallel configuration can be restored, when an exter- 
nal magnetic field, which compensates the exchange field 
created by the ferromagnetic leads, is applied] 10 ' 11 This 
behavior was confirmed in a couple of recent experi- 
ments J^— 

The situation becomes more complex and physically 
richer when the dot is exchange-coupled to an additional 
reservoir J£ Such a model captures the essential physics of 
the so-called two-channel Kondo effect.— ~— In the two- 
channel Kondo problem, two separate electron reservoirs 
(channels) compete with each other to screen the im- 



purity's spin. If the coupling to one of them is larger 
than to the other one, a usual single-channel Kondo 
state (spin singlet) is formed between the dot and more 
strongly coupled reservoir. This results in two compet- 
ing Kondo ground states of the system, depending on the 
ratio of coupling strengths to the first and second conduc- 
tion channels. Interestingly, these two Kondo states are 
separated by a quantum critical point, where both cou- 
plings are equal and an exotic two-channel Kondo state 
is formed, which cannot be described within the Lan- 
dau Fermi-liquid theory. Very recently, the two-channel 
Kondo effect has been explored experimentally in quan- 
tum dots^ The experimental setup consisted of a small 
quantum dot coupled to external leads and to a large 
Coulomb-blockaded island. While the electrons could 
tunnel between the dot and the leads, only virtual tun- 
neling processes between the dot and the island were al- 
lowed, resulting in an exchange coupling. By tuning the 
exchange coupling, it was possible to study the quantum 
phase transition between the two ground states of the 
system and analyze transport behavior in the non-Fermi 
liquid regime^ Theoretically, such a two-channel setup 
can be modelled for example by a quantum dot which 
is tunnel-coupled to external leads and exchange-coupled 
to another electron reservoir 

As discussed above, both the Kondo effect in a quan- 
tum dot coupled to ferromagnetic leads and the two- 
channel Kondo phenomenon were already extensively 
studied. However, the interplay of leads' ferromagnetism 
and two-channel Kondo effect remains to a large extent 
unexplored. Therefore, in this paper we address the 
two-channel Kondo problem in the presence of ferromag- 
netism. In particular, we consider an Anderson quantum 
dot coupled to ferromagnetic leads and exchange-coupled 
to a nonmagnetic electron reservoir. Using the numerical 
renormalization group (NRG) method, we analyze the in- 
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FIG. 1: (color online) The schematic of a quantum dot (QD) 
tunnel-coupled to external ferromagnetic leads and exchange- 
coupled to a nonmagnetic electron reservoir. The spin- 
dependent coupling to the left (right) lead is described by 
Tlo- (Trct), while J denotes the exchange coupling constant. 
The magnetizations of the leads can form either parallel or 
antiparallel magnetic configuration, as indicated. 

tcrplay between the effects due to ferromagnctism of the 
leads and exchange coupling to the additional nonmag- 
netic reservoir. Depending on the strength of the tunnel 
coupling t and exchange coupling J, the dot's spin can 
be screened either by electrons in the ferromagnetic leads 
or by electrons in the nonmagnetic reservoir. By analyz- 
ing the equilibrium spectral functions and the conduc- 
tance through the dot, we show that in the antiparallel 
magnetic configuration, the system undergoes a quan- 
tum phase transition with increasing exchange coupling 
J, where the conductance drops from the maximum value 
to zero. For a certain critical value of J, J^ p , both elec- 
tron channels try to screen the dot's spin and the con- 
ductance approaches a half of the quantum conductance. 
In the parallel configuration, on the other hand, the con- 
ductance is generally suppressed, irrespective of the ex- 
change constant J, due to effective spin splitting of the 
dot level caused by the exchange field coming from fer- 
romagnetic leadsi&£ We show that the Kondo effect can 
be restored by applying a properly tuned external mag- 
netic field B for J below the corresponding critical point, 
J < . Furthermore, the quantum critical regime can 
also be recovered, which however requires a fine-tuning 
in the parameter space of J and B. 

The paper is organized as follows. In section II we 
present the model as well as briefly describe the NRG 
method together with some details of calculations. In 
turn, in section III we present numerical results for sym- 
metric and asymmetric Anderson models in both parallel 
and antiparallel magnetic configurations of the system. 
Finally, we conclude in section IV. 

II. THEORETICAL DESCRIPTION 
A. Model 

The considered system consists of a single-level quan- 
tum dot tunnel-coupled to left and right ferromagnetic 
leads and exchange-coupled to a nonmagnetic reservoir, 



see Fig. [I] It is assumed that the external leads are made 
of the same ferromagnetic material and their magnetiza- 
tions are collinear, so that the system can be either in 
the parallel or antiparallel magnetic configuration. The 
total Hamiltonian is given by 

H = HpM + ^NM + Hqv + H tu n + -ffcxch- (1) 

Here, Hfm describes the ferromagnetic leads, iJpM = 
J2rka e rkaC rka c r ka, where c\ kry is the electron creation 
operator with wave number fc, spin a in the left [r = L) or 
right (r = R) lead, and e r ka is the corresponding energy. 
The second part, -J/nm, corresponds to a nonmagnetic 
electron reservoir and is given by, J?nm = J2ka £ k a \a a kcn 
with a\ a being the respective creation operator and £k is 
the single-particle energy. The quantum dot is described 
by the Anderson Hamiltonian, 

i? QD = J2 £dd ° d ° + ^4^4^ + BS * - ( 2 ) 

(7 

where d\ creates a spin-cr electron, Ed denotes the en- 
ergy of an electron in the dot, and U describes the 
Coulomb correlations between two electrons occupying 
the dot. The last term corresponds to external magnetic 
field B applied along the zth direction {gps = 1) an d 
5*2 = 5(4 ( ^1" — 4^)' tunnel Hamiltonian is given 
by 

#tun = ^ tra (dtcrka + c lka d <r) J ( 3 ) 
rka 

where t ri7 describes the spin-dependent hopping matrix 
elements between the dot and ferromagnetic leads. The 
coupling to magnetic leads can be described by Y ra = 
■np r \t riJ \ 2 , where p r = p is the density of states in the 
lead r. We have thus shifted the whole spin-dependence 
into the coupling constants and assumed a flat band of 
width 2D^2. where D = 1 is set as the energy unit, if 
not stated otherwise. 

By means of a unitary transformation in the left-right 
basis ; 31 ! 32 one can map the problem of tunneling through 
quantum dot coupled to the left and right leads into a 
problem where the dot is effectively coupled to a single 
lead with a new coupling constant, T a = Tlct + Tro-. 
This can be done by introducing the following symmetric 
operators, Ofc CT = t ha c Lkl ? + t Rcr CR fc a -, with dimensionless 
coefficients t r<7 = t ra / y/ 1 2 ct + t^ a . Then, the tunneling 
Hamiltonian can be written as 

k(7 V 

One can see that now the dot is tunnel-coupled to only 
one effective electron reservoir, £/fm = J2ka £ ka- a k a a k<r, 
with new spin-dependent coupling constant T a . The 
other parts of the system Hamiltonian, Eq. (TIJ , are not 
affected by this transformation. To parameterize the 
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spin-dcpcndcnt couplings we also introduce the spin po- 
larization of ferromagnetic leads, p = (T^ — T±)/(T^+r±). 
The couplings can be then written in a compact form as, 
r m) = (1 ±p)T, where T = (r T + Assuming 
symmetric coupling strength of the dot to the leads, the 
resultant coupling in the antiparallel configuration is the 
same for the spin-up and spin-down electrons, = V. 
On the other hand, in the parallel configuration, the 
couplings are then different for the two spin directions, 
I^rn = (l±p)r, which effectively leads to spin splitting 
of the dot level and, when this is the case, the Kondo res- 
onance may become suppressed because of broken spin 
degeneracy^ 

Finally, the exchange Hamiltonian describing the cou- 
pling between the dot and the second (nonmagnetic) 
reservoir is given by 



<a.ka 



(5) 



where S — \ X^o-o-' d\,a aa id a i is the spin in the dot, J de- 
notes the exchange coupling constant and a is a vector of 
Pauli spin matrices. We note that in addition to the ex- 
change scattering of electrons [described by Eq. (5)] there 
could be also potential scattering. However, in this work 
we are mainly interested in the low energy physics, where 
the Kondo effect emerges, so the potential scattering may 
be neglected, as it docs not lead to any Kondo-type cor- 
relations. 



both Abelian and non-Abelian symmetries! 35 ' 36 In calcu- 
lations we have thus used the U(l) symmetry for the zth 
component of the total spin, the £7(1) symmetry for the 
charge in the first channel, and the SU(2) symmetry for 
the charge in the second channel. Furthermore, in calcu- 
lations we have taken the discretization parameter A = 2 
and kept 3000 states at each iteration step. 

Using the NRG we can calculate the spectral function 
of the dot, A a {oS) = -±ImG% a {w), where Gf a (w) de- 
notes the Fourier transform of the dot retarded Green's 
function, G% a (t) = -ie(t){{d a (t), 4(0)}). On the other 
hand, the spectral function can be directly related to the 
spin-resolved linear conductance G a by the following for- 
mula 



e 2 4r LCT r RCT 



duntA a {uj) 



dui 



(6) 



where /(w) is the Fermi distribution function and the to- 
tal conductance is given by, G — Gf + G\. At zero tem- 
perature, the spin dependent conductance for the paral- 
lel configuration is given by, G^m = ± P^FA^y 
while for the antiparallel configuration one gets, G^^ = 

5^(1 - p^irTA^ = G AP /2, where A^ /AP is the zero- 
temperature spectral function of the d-level operator in 
respective magnetic configuration, taken at w = 0. 



B. Method 

To analyze the equilibrium transport properties of the 
considered system, we employ the numerical renormal- 
ization group method^ 3 - - nonperturbative, very powerful 
and essentially exact numerical method to address quan- 
tum impurity problems.— The NRG consists in a loga- 
rithmic discretization of the conduction band and map- 
ping of the system onto a semi-infinite chain with the 
impurity (quantum dot) sitting at the end of the chain. 
By diagonalizing the Hamiltonian at consecutive sites of 
the chain and storing the eigenvalues and eigenvectors 
of the system, one can calculate the static and dynamic 
quantities of the system. In the case of model consid- 
ered in this paper, the Hamiltonian is mapped onto two 
semi-infinite chains, where the first chain corresponds to 
ferromagnetic leads tunnel-coupled to the dot, while the 
second one to nonmagnetic reservoir exchange-coupled 
to the dot. Because such two-channel calculations are 
usually very demanding numerically, it is crucial to ex- 
ploit as many symmetries of the system's Hamiltonian 
as possible. Especially, using the ST/ (2) symmetry de- 
creases the size of Hilbert space and thus increases con- 
siderably the accuracy of calculations. In particular, to 
efficiently perform the analysis, we have used the flexible 
density-matrix numerical renormalization group (DM- 
NRG) code, which can tackle with arbitrary number of 



III. NUMERICAL RESULTS 

In the following we present numerical results on the 
equilibrium spectral function and linear conductance, 
when the quantum dot is in the Kondo regime. We will 
distinguish between two different situations; symmetric 
(e<j = —U/2) and asymmetric (£<j ^ — U/2) Anderson 
models. The origin of such a distinction stems from the 
way in which ferromagnetic leads act on the quantum dot. 
More specifically, in the asymmetric Anderson model fcr- 
romagnetism of the leads gives rise to a spin splitting of 
the Kondo resonance in the parallel configuration, while 
in the symmetric model no such a splitting appears (as- 
suming that the dot is coupled with the same strength to 
the left and right leads) £iLOL In other words, an effective 
exchange field, due to coupling to magnetic leads, acts 
on the dot in the former case, while such a field vanishes 
in the latter case. The effective field is directly related 
to the difference in the coupling strengths of the dot and 
ferromagnetic leads for the two spin orientations. Since 
the coupling in the spin-up channel is larger than that 
in the spin-down one, energy of the spin-up (spin-down) 
electron in the dot decreases (increases) by Asa/2. Con- 
sequently, the spin-dependent coupling acts as an effec- 
tive magnetic field, leading to spin-splitting Asd of the 
dot level££±2 
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FIG. 2: (color online) The spectral function of the d-level 
operator in the antiparallel (a) and parallel (b,c) magnetic 
configurations for the symmetric Anderson model and for dif- 
ferent values of the exchange coupling J. The parameters 
are: e d = -0.05, U = 0.1, T = 0.0077, p = 0.4, and T = 0. 
The Kondo temperature is defined as a half-width of the spec- 
tral function for J = and p = 0, Tk = 2.5 x 10 -4 , while 
Aa = J^ct Acr(La — 0) for J = and p = 0. All the parameters 
are given in the units of D = 1. 



A. Symmetric Anderson model 



For the symmetric Anderson model we assume the fol- 
lowing parameters (in the units of D), s c \ = —0.05 and 
U = 0.1. The zero-temperature spin-dependent spectral 
function A a , normalized to Aq, with Aq = A a {uj = 0) 
taken for J = and p = 0, is shown in Fig. [5]for both an- 
tiparallel (a) and parallel (b,c) magnetic configurations 
(note the logarithmic energy scale) , and for indicated val- 
ues of the exchange coupling parameter J. The spectral 
function is plotted as a function of uj/Tk, where Tk is the 




FIG. 3: The linear conductance as a function of exchange cou- 
pling constant J for the antiparallel magnetic configuration. 
The conductance was determined from the spectral function 
shown in Fig. [21 a). The parameters are the same as in Fig. [5] 
and J is in units of D — 1. 



Kondo temperature defined as a half-width of the d-level 
spectral function for J = and p = 0, Tk = 2.5 x 10~ 4 . 
It can be seen that in the antiparallel configuration the 
spectral function is independent of the spin orientation 
[Fig. [H[a)], A^ p = A^ p , while it depends on electron 
spin in the parallel magnetic configuration, A P ^ A P , 
see Fig. &[b) and (c). Note, that for symmetric Ander- 
son model the spectral function is symmetric with respect 
to uj = 0, therefore here it is shown only for positive en- 
ergies, i.e. for energies above the Fermi level. Moreover, 
note also that the spectral functions are normalized to 
that for the corresponding paramagnetic limit (p = 
and J = 0), so A^(u = 0) + A±(uj = 0) ^ A in the 
parallel configuration. 

Let us consider first the situation with vanishing ex- 
change coupling of the dot to the nonmagnetic reservoir, 
J = 0. For uj < Tk, a Kondo peak develops in the dot 
spectral function due to screening of the dot's spin by 
conduction electrons of ferromagnetic leads, which leads 
to the formation of a non-local spin singlet. The height 
of the Kondo peak is independent of spin in the antipar- 
allel configuration and depends on spin in the parallel 
one. Apart from this, a Hubbard peak corresponding to 
£d + U is visible in the spectral function shown in Fig. [21 
This behavior of the Kondo phenomenon in the presence 
of ferromagnetic leads is in agreement with that found 
by other methods, for instance by the equation of mo- 
tion for the Green functions^ and also by the real-time 
diagrammatic technique^ 

The situation changes when the electron in the dot is 
additionally exchange coupled to the nonmagnetic reser- 
voir. When the coupling is antiferromagnetic and the 
coupling parameter J increases, the width of the Kondo 
peak becomes gradually narrower and narrower. The 
hight of the peak, however, remains unchanged, as can 
be clearly seen in Fig. [5] for some small values of the ex- 
change coupling constant. In order to see this behavior 
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also for larger J, but still smaller than a critical value, 
one should plot the spectral function for 

P(AP) 

lower energies. For J < J c , the system is in the spin 
singlet ground state formed by the quantum dot spin and 
electrons in the ferromagnetic leads, which gives rise to 
the Kondo resonance in the spectral function. However, 
when J > Jc , the coupling to nonmagnetic reser- 
voir becomes larger than the coupling to ferromagnetic 
leads and the dot's spin becomes screened by electrons of 
the nonmagnetic reservoir. Now the Kondo peak in the 
spectral function disappears for both magnetic configu- 
rations of the system, see Fig. [2] When the two couplings 
are equal, i.e. for J = J^ AP '', the system is in an ex- 
otic state where the two channels try to screen the dot's 
spin. The spectral function at u = is then equal to 

1 A P /AP 

a half of its value corresponding to J = 0, 2 At j=o, 
for both magnetic configurations, see Fig. [2j This be- 
havior reveals a quantum phase transition with increas- 
ing strength of the exchange coupling. The origin of the 
phase transition follows from the interplay of the tunnel 
coupling to the ferromagnetic electrodes and exchange 
coupling to the nonmagnetic reservoir. More specifically, 
the quantum phase transition occurs at the boundary be- 
tween two different singlet ground states, involving the 
dot's spin and conduction electrons of the leads or side- 
coupled reservoir. The behavior of transport characteris- 
tics around this critical point in the case of nonmagnetic 
system was discussed in Ref. [28| . It was shown that the 
zero-temperature conductance depends step-like on the 
difference A between the tunnel and exchange couplings, 
and becomes equal to a half of its maximum value at 
the critical point, i.e. when A = 0. The discontinu- 
ity of the linear conductance with respect to A reflects 
the quantum phase transition in the parameter space of 
tunnel coupling t and exchange coupling J. Since the 
conductance is determined by the corresponding spectral 
functions at w = 0, quantum critical behavior is also 
reflected in the J-dependence of the spectral function. 
We also note that at finite temperature the transition is 
smeared as \/T/Tk and turns rather into a crossover^ 

Using the Schricffcr- Wolff transformation^ one could 
try to estimate the critical value of J. For the symmetric 
Anderson model and for antiparallcl configuration one 
gets^ J c s - W = r(7r / 9)- 1 [/|£ d |- 1 (e d + [/)^ 1 « 0.196. From 
the numerical data, however, one finds J C AP « 0.22878 
for the antiparallcl and J P « 0.22138 for the parallel 
configurations, see Fig. [2] The difference between the 
value obtained using the Schricffcr- Wolff transformation 
and the numerical value may result for example from 
the fact that the transformation is based on perturbation 
expansion, and takes into account only the second-order 
tunneling processes. 

The quantum critical behavior can be also seen in the 
dependence of the linear conductance on the coupling 
constant J, which is shown in Fig. [3] for the antiparallcl 
magnetic configuration. For J < J AP , the conductance 
is G AP = (1 — p 2 )2e 2 /h and drops to zero when J > 




P 



FIG. 4: (color online) The dependence of the critical exchange 
coupling (in units of D = 1) on the spin polarization of the 
leads p in the case of symmetric Anderson model and parallel 
magnetic configuration for three different values of the tunnel 
coupling r, as indicated in the figure. The other parameters 
are the same as in Fig. [2] 



. On the other hand, at the quantum critical point 
J = J AP , the linear conductance is equal to half of its 
value for J = 0, i.e. G AP = (1 — p 2 )e 2 /h. Consequently, 
the dependence of the conductance on the exchange cou- 
pling J can be expressed as G AP = 0(A)(1 — p 2 )2e 2 /h, 
where A = J AP — J and Q(x) is the step function. The 
dependence of G on J for the parallel configuration is 
qualitatively similar to that in the antiparallel configura- 
tion, therefore it is not shown here. 

In Fig. H]the dependence of the critical exchange cou- 
pling J P on the spin polarization of the leads p in the case 
of symmetric Anderson model and parallel magnetic con- 
figuration is shown for three different values of the tunnel 
coupling T. First of all, the critical coupling J P decreases 
with decreasing the coupling strength T. Moreover, J P 
also decreases with increasing the spin polarization of the 
leads. For p — > 1, J P tends to zero, as only spins of one 
orientation are coupled to the leads and the Kondo effect 
becomes suppressed. This behavior of the critical param- 
eter J P is consistent with the dependence of the Kondo 
temperature in a quantum dot coupled to ferromagnetic 
leads on the coupling strength T and spin polarization 



B. Asymmetric Anderson model 

Let us now consider the case of asymmetric Anderson 
model, |ea| 7^ U/2. For numerical calculations we as- 
sume £<j = —0.05 and U = 0.2. The spectral function 
in the antiparallel magnetic configuration as a function 
of lu/T k , where T K = 3.4 x 10~ 5 , is shown in Fig. [5] for 
indicated values of the exchange coupling parameter J. 
The inset shows the behavior of the spectral function as- 
sociated with the Kondo peak. The general features of 
the spectral function are similar to those of the corre- 
sponding spectral function in the case of symmetric An- 
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FIG. 5: (color online) The spectral function of the d-level 
operator in the antiparallel magnetic configuration for the 
asymmetric Anderson model, £d = —0.05, U — 0.2, and for 
indicated values of the exchange coupling J. The Kondo tem- 
perature for assumed parameters (and for J = and p = 0) 
is Tk = 3.4 x 10 -5 . The other parameters are the same as in 
Fig-H 



derson model discussed above, see Fig. [2] This is because 
in the antiparallel configuration the resultant coupling to 
ferromagnetic leads does not depend on spin and the sys- 
tem effectively behaves as a nonmagnetic one. As before, 
one observes a quantum phase transition at J = , 
where now J^ p « 0.1858. The only difference is that for 
|e<j| 7^ U/2 the spectral function displays an asymmetric 
behavior with respect to U) = 0, see the inset in Fig. [SJ 

The situation, however, changes significantly when the 
magnetizations of the leads switch to the parallel config- 
uration. The corresponding spectral function for spin-t 
and spin- J, is shown in Fig. [5] Note, that now the spectral 
function is shown for both positive and negative energies. 
As before, let us consider first the case of J — 0. Due to 
an effective exchange field originating from the presence 
of ferromagnetic electrodes, the spin degeneracy of the 
dot level is lifted. At zero temperature, the magnitude of 
the splitting due to exchange field, Aed, can be estimated 
from the formula£ii£ 



Ae ri = 



2pT 



7T \s d + U\ 



(7) 



For the assumed parameters one then finds, Aed « 
2.15 x 10 -3 . The exchange field leads generally to the 
suppression of the Kondo peak, however the reminiscent 
of the Kondo effect are still visible as relatively small 
peaks in the d-level spectral function. The position of 
these peaks is shifted away from the Fermi level - to 
positive energies for spin- J, and to negative energies for 
spin-t- In fact, the peaks occur for energies comparable 
to the magnitude of the exchange field. As can be seen 
in Fig. |B1 they develop at w/Tr- ps 10 2 for spin- 4- and 
at lu/Tk ~ — 10 2 for spin-j" components of the spectral 
function. The other peaks in the spectral functions cor- 




-10 4 -10 3 -10 2 -10 1 10 1 10 2 10 3 10 4 
w/T K 

FIG. 6: (color online) The spectral function of the d-level 
operator in the parallel magnetic configuration for the asym- 
metric Anderson model. The other parameters are the same 
as in Fig. [5] 



respond to the dot level £d and its Coulomb counterpart 
e d + U. 

When the coupling parameter J increases, the weak 
Kondo resonances in the spectral function gradually dis- 
appear for both spin orientations. The physics behind 
this disappearance remains similar to that described 
above, i.e. screening of the dot's spin by the nonmag- 
netic reservoir exchange-coupled to the dot. Interest- 
ingly, there is no quantum phase transition in the case of 
parallel magnetic configuration shown in Fig. [5] 

Let us now assume that there is an external magnetic 
field B applied to the dot along the zth direction. In 
the case of antiparallel configuration, the magnetic field 
destroys both the Kondo resonance and quantum phase 
transition with changing J. However, when the leads are 
aligned in parallel, the Kondo effect is already suppressed 
by the effective exchange field coming from ferromagnetic 
electrodes, and one may consider the possibility of restor- 
ing the Kondo peak by applying an external magnetic 
field which compensates the effects due to exchange field. 
In Fig. [7] we show the spectral functions for the parallel 
magnetic configuration in the case of an asymmetric An- 
derson model, calculated for three different values of the 
exchange constant J in the presence of external magnetic 
field B. The insets display behavior of the spectral func- 
tion associated with the Kondo peaks. When J = 0, see 
Fig.JTJa) and (b), the full Kondo peak at the Fermi level 
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FIG. 7: (color online) The spectral function of the d-level operator in the parallel magnetic configuration for the asymmetric 
Anderson model in the presence of external magnetic field B applied along the zth direction for different values of exchange 
coupling constant J = (a,b), J = 0.17 (c,d) and J = 0.17875 (e,f). The other parameters are the same as in Fig. [5] 



in the spectral density can be restored for both spin orien- 
tations by properly tuned external magnetic field, which 
happens for B = B c = 0.0020545, where B c (in the units 
of D) denotes the compensating field. This is in agree- 
ment with the result obtained earlier Similar behavior 
also appears for larger positive J, e.g. for J = 0.17 shown 
in Fig.[7Jc) and (d). Now, the Kondo resonance becomes 
restored when the compensating field is B c = 0.001928. 
Note that the magnitude of magnetic field necessary for 
full restoration of the Kondo effect slightly decreases as J 
increases. The question which arises now is whether such 
a restoration by magnetic field is also possible for larger 
values of J. By fine-tuning in the parameter space of J 
and B, we have found that this is the case for J below 
a certain critical value J < = 0.17875. Here J<f de- 
notes the critical value of J in the parallel configuration 
and in the presence of the compensating magnetic field. 
From numerical results (not shown here), follows that for 
J > J^, the magnetic field can only partially restore the 
Kondo effect, leading to small side peaks in the spectral 
function, while the full Kondo peak at tu = cannot be 
restored. One may now expect that for J = jf , the mag- 
netic field should also restore the quantum critical state. 
Indeed, by fine-tuning in the parameter space we have 
found that the quantum critical state can be recovered 
for B c = 0.00191621. This situation is shown explicitly 
in Fig. [T^e) and (f). Thus, we have shown that in the 
parallel configuration the properly-tuned magnetic field 
can restore both the full Kondo effect for J < as well 
as the quantum critical state for J = . 

By comparing numerical curves presented in Fig. one 



0.00206 




FIG. 8: The dependence of the compensating magnetic field 
B c on the exchange coupling constant J in the case of the 
parallel configuration and asymmetric Anderson model. The 
other parameters are the same as in Fig. [5] J and B c are in 
units of D = 1. 



can note that the compensating field B c decreases with 
increasing the exchange coupling J. This is explicitly 
shown in Fig. [8l where we have calculated the depen- 
dence of B c on the exchange coupling J. For J < , 
the Kondo resonance can be fully restored by applying 
compensating field B c . On the other hand, when J > , 
the magnetic field cannot compensate the exchange field, 
so the notion of compensating field becomes meaningless. 
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IV. CONCLUSIONS 

In this paper we have considered spectral and trans- 
port properties of a single-level quantum dot connected 
to external ferromagnetic leads and exchange-coupled to 
a nonmagnetic reservoir. Using the numerical renor- 
malization group method we have calculated the zero- 
temperature d-level spectral function and the conduc- 
tance through the dot. We have shown that in the an- 
tiparallel configuration, depending on the strength of the 
exchange interaction J, the Kondo singlet ground state 
can form, in which the conduction electrons either in the 
ferromagnetic leads or in the nonmagnetic reservoir are 
involved. In the former case, the conductance is maxi- 
mum, whereas in the latter case the conductance becomes 
fully suppressed. For a certain critical value of J, J AP , 
both electron channels try to screen the dot's spin and 
the conductance is equal to a half of its maximum value. 
The boundary between the two ground states is a quan- 
tum phase transition. 

In the parallel magnetic configuration, on the other 
hand, the Kondo effect is generally destroyed due to an 
effective spin splitting of the dot level caused by the pres- 
ence of ferromagnetic leads. However, there are still small 
side peaks - reminiscent of the Kondo effect - which oc- 
cur on both sides of the Fermi level for energies of the 
order of effective exchange field. Nevertheless, with in- 



creasing the exchange constant J, these peaks become 
suppressed. 

We have also considered the influence of an external 
magnetic field on the d-level spectral function and shown 
that in the parallel configuration the Kondo effect can be 
restored by applying appropriately tuned compensating 
magnetic field for J < J P , where J P is the critical value 
of J in the compensating magnetic field. If, however, 
J > J P , the full Kondo effect cannot be restored by a 
magnetic field. In addition, we have found that the quan- 
tum critical behavior, which is suppressed in the parallel 
configuration, can also be recovered by tuning the exter- 
nal magnetic field. 
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